Discovery of Relativistic Outflow in the Seyfert Galaxy Ark 564 
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ABSTRACT 

We present Chandra high energy transmission grating spectra of the narrow- 
line Seyfert- 1 galaxy Ark 564. The spectrum shows numerous absorption lines 
which are well mod eled with low veloc ity outflow components usually observed 



in Seyfert galaxies (IGupta et al.ll2013l ). There are, however, some residual ab- 
sorption lines which are not accounted for by low-velocity outflows. Here we 
present identifications of the strongest lines as Ka transitions of O Vll(two lines) 
and Ovi at outflow velocities of ~ 0.1c. These lines are detected at 6.9a, 6.2a, 
and 4.7cr respectively and cannot be due to chance statistical fluctuations. Pho- 
toionization models with ultra-high velocity components improves the spectral 
fit significantly, providing further support for the presence of relativistic outflow 
in this source. Without knowing the location of the absorber, its mass and en- 
ergy outflow rates cannot be well constrained; we find E (out flow) /L bo ranging 
from < 0.001% to 60% using different assumptions. This is the first time that 
absorption lines with ultra-high velocities are unambiguously detected in the soft 
X-ray band. The presence of outflows with relativistic velocities in AGNs with 
Seyfert-type luminosities is hard to understand and provides valuable constraints 
to models of AGN outflows. 
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Introduction 



There have been reports of ultra-high velocity outflows in the X-ray spectra of radio- 
quiet AGNs and quasars. These outflows are identified through blue-shifted Fexxv and/or 
Fexxvi absorption lines (at rest-frame energies of 7-10 keV) from highly ionized gas {log^ = 
3 — 6 erg s _1 ), with col umn densities as larg e a s Nfj = 10 22 — 10 24 cm " 2 , and with relativistic 
velocities of 0.1c-0.3c ( Pounds et al. 2003a fbL ITombesi et al. 2010 , and references therein). 
The mass outflow rate of these high velocity outflows were comparable to the accretion rate 



and t heir kinetic energy was a significant fraction of the bolometric luminosity (ITombesi et al. 



20121 ). These outflows can provide the effective feedback required by theoretical models of 
galaxy formation to solve a number of astrophysical problems ranging from cluster cooling 
flows to structures of galaxies. 

Although high velocity outflows are detected in a large number of sources, in all the 
cases identification are based on Fe K-shell transitions which fall in a region of the spectrum 
where instrumental resolution is much lower than in the soft X-ray band. The significance of 
the absorption line detections is often questioned and with only a few lines observed, accurate 
parametrization of the photoionized plasma becomes difficult. While the observed transitions 
are from highly ionized gas (log£ = 3 — 9 ers s -1 ), photoionization models predict that in such 
a plasma in addition to highly ionized iron, lighter elements should a lso be present, strongest 
of which are S, Si and O (e.g.. lSim et al.ll2010l . iKrongold et al.ll2003l ). Absorption lines from 
these highly ionized elements lie in the soft X-ray band, so warm absorber signatures of 
relativistic outflows at soft X-rays should be present, but were never unambiguously observed 
until now (see §5.3 for existing evidence). For complete understanding of the properties of 
high velocity outflows it is necessary to detect transitions of other ions such as of O, Ne, Si 
and S as predicted by models. 

In this paper, we report the serendipitous discovery of high velocity outflows in the soft 
X-ray band in Ark 564, identified during our detailed analysis of the Chandra archival data 
of this source ( jGupta et al.ll2013l hereafter Paper I). 



Primary goal of Paper I was to self consistently analyze and model the grating spectra 
of typical (V ou t = 100 — 1000 km s^ 1 ) warm absorbers in Ark 564. Here we present the 
discovery of warm absorbers (WA) with ultra-high velocities. 



2. Data and Spectral Analysis 



Ark 564 is a bright, nearby (z = 0.024684), narr ow-line Seyfert 1 (NLS1) galaxy, with 



luminosity of L2- 



10 keV 



(2.4 - 2.8) x 10 43 ergs s" 1 (ITurner et aI.ll200lL iMatsumoto et al. 
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2004j . Paper I). In Paper-I we have discussed the Chandra observation and data reduction of 
this source which we briefly summarize here. Ark 564 was observed with the Chandra High 
Energy Transmission Grating Spectrometer (HETGS) in 2008 in three separate exposures 
totaling 250 ks. We followed the standard procedure to extract spectra. We used the software 
package CIAO (Version 4.3) and calibration database CALDB (Version 4.4.2) developed by 
the Chandra X-ray center. We co-added the negative and positive first-order spectra and 
built the effective area files (ARFs) for each observation using the fullgarf CIAO script. The 
HETG-MEG spectra were analyzed using the CIAO fitting package Sherpa. 

We co-added the spectra obtained with HETG-MEG and averaged the associated ARFs. 
We fitted the continuum with absorbed power law plus black body component and modeled 
all the statistically significant local absorption (z = 0) features with Gaussian components. 
Further we used the photoionization m odel fitting code PHotoionized Absorption Spectral 



Engine (PHASE; iKrongold et al.l 120031 ) . to fit the typical warm absorber features. The best 
fit model of intrinsic absorption requires two warm absorbers with two different ionization 
states (logU = 0.38 ± 0.02 and loglJ = —1.3 ± 0.13), both with moderate outflow velocities 
(~ 100 km s" 1 ) and relatively low line of sight column densities (~ Nh = 10 20 cm" 2 ). For 
detailed spectral analysis and best fit parameters, we refer readers to Paper-I. 

Though most of the intrinsic absorption features in the source spectrum are well fitted 
with two warm absorbers, there are residual absorption line-like features in the spectral 
regions of 19.0 — 21.0 A, 17.0 — 17.5 A and 13.5 — 13.8 A each with individual significance 
of 1.7 — 6.9 sigma (Fig. 1). To check for the consistency of these residuals and to confirm 
that these are not the artifacts of co-adding spectra, we inspect individual Ark 564 HETGS- 
MEG spectra from the 2008 observations. After fitting the data with continuum, Galactic 
and two-phase WAs model (Model-A) as noted above, the same residuals are found in the 
same spectral region in all the three observations (Fig. 2). In the following sections, we 
discuss the identifications, statistical significance and possible origin of these absorption 
lines. 



3. Discovery of the high velocity outflows 

3.0.1. Identification of the Absorption Lines 

The strongest residual features in the co-added spectrum are present at 19.805±0.005 A, 
19.845 ±0.001 A and 20.250 ±0.011 A (observed frame, Fig. 1). Errors refer to 1 a confidence 
level throughout the paper, unless noted otherwise. First we try to identify the features 
at 19.805 ± 0.005 A and 19.845 ± 0.001 A. There is no known instrumental feature near 
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these energies (Chandra Proposers' Observatory Guide, or POG). The z=0 lines are already 
included in the model; there are no permitted lines with oscillator strength > 0.0001 at 
wavelength of either 19.805, 19.845 or 20.250 A. At the observed energies, there would be no 
intervening system with Zwhtm < zjrk rr 4 from the warm- hot intergalactic medium (WHIM; 



Mathur et al.ll2003l . iNicastro et al. II2005I ). Therefore, we assume that these absorption lines 
are intrinsic to the source. We identify these lines based on a combination of chemical 
abundance and line strength (the oscillator strength / > 0.1). Given the detected wavelength 
and assuming a very broad range of inflow/outflow velocities —60, 000 to 60, 000 km s -1 the 
likely candidates are Ovil K/3 at \ rest = 18.62A, Ovill Ka at \ rest = 18.96A, Caxvil k/3 
at Kest = 19.56A, Caxvi k/3 at \ rest = 20.617A, Arxiv k/3 at X rest = 21.15A, Caxvi Ka 
at Xrest = 21.45A and Ovil ka at X re st = 21.602A. Considering that argon and calcium 
are orders of magnitude less abundant than oxygen, and because O VII Ka and O vm Ka 
are by far the strongest possible lines, the most likely candidates are O vm Ka with inflow 
velocities of 0.019c and 0.021c, or O VII Ka with outflow velocities of 0.105c and 0.103c. 

To distinguish between the two possibilities of inflow and outflow, we search for possible 
associations of other lines such as Ovill kj3, Ovil k/3, Ovi and/or Ov . We do not find 
any possible association for inflows. The absorption feature at 20.25 ± O.OlA corresponds to 
Ovi ka (X rest = 22.037A) at the outflow velocity of 0.103c. We also found an absorption 
feature at 17.08 ±0.01 A with 1.7a significance, corresponding to O VII k/3 line at the outflow 
velocity of 0.105c. The detection of Ovi Ka and O VII K/3 lines at the same velocity favors 
the outflow scenario. 



3.0.2. Statistical Significance of Absorption Lines 

To investigate whether the apparent absorption lines could be due to statistical fluctu- 
ations, we calculate the probability of detection of individual lines due to random statistical 
fluctuation. First we fit the lines with negative gaussians of fixed width 0.001A, folding 
through the detector response and leaving other parameters (centroid and amplitude) free 
to vary. The addition of three gaussian lines at 19.805 ± 0.005 A, 19.845 ± 0.001 A and 
20.250 ± 0.011 A to our previous model (Model A) improves the fit statistic by Ax 2 = 
120, Ad.o.f. = 6, an improvement at more than 99.99% confidence by the F-test (Bevington 
and Robinson 1992). We measured the equivalent width (EW) of lines at 19.805 A, 19.845 A 
and 20.250A of 15.6 ±2.5 mA, 16.5 ± 2.4 mA and 16.1 ±3.4 mA respectively. Errors are at 
1 sigma confidence level and are calculated using the "projection" command in Sherpa. Thus 
the three absorption lines are detected with 6.2o~, 6.9cr and 4.7<r significance respectively. 
Further using the Gaussian probability distribution, we looked for the probability of detec- 
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tion of these lines by chance. For the lines detected with 6.2a, 6.9a and 4.7a significance, 
the probability of false detection is 2.8 x 10~ 10 , 2.6 x 10~ 12 and 1.3 x 10~ 6 respectively. There 
are, however, 4801 wavelength bins in our spectrum (these are all the wavelength bins in the 
wavelength range of 1 — 25 A, beyond which data is of poor quality. The rest of the HETG 
spectrum was never used in any of the analysis, even in Paper I). Therefore the probability 
of finding absorption lines at the observed significance anywhere in the spectrum due to 
random statistical fluctuations is 0.0001%, 0.000001% and 0.6% respectively. 

For the outflow system with velocity of 0.105c, we also detected the O VII Kf3 line with 
chance probability of 0.04. Thus the combined chance detection probability of detecting 
both Ovil Ka and Ovn K(3 lines is 4 x 10 _6 %. Similarly for the other system (at 0.103c) 
the combined chance detection probability of detecting both O VII Ka and O VI Ka is 4 x 
10 _8 %. Thus we conclude that the detected absorption lines are not due to random statistical 
fluctuations, but are signatures of outflows. 

3.0.3. Other Absorption lines at blueshift o/~ 0.1c 

After confirming the high detection significance of O VII and O VI absorption lines iden- 
tified with high velocity outflows, we searched the entire spectrum for other ionic transitions 
at the similar blueshifts. We detected two other absorption features at 17.35 ± 0.01 and 
13.63 ± 0.01 A with EWs of 11.8 ± 2.7 and 5.0 ± 1.7 ml (Fig. 3, Table 1). The feature at 
13.63A is also detected in the HEG spectrum. The tentative identification of these features 
are O vm Ka and Fexvil Ka at outflow velocities (0.110 ± 0.001)c and (0.114 ± 0.001)c 
respectively. The outflow velocities of these systems are close to each other and to the above 
reported systems of high velocity outflow. 

4. Photoionization Modelling of High Velocity Outflows 

To determine the physical properties of the absorber responsible for producing highly 
blueshifted absorption features and to check for the physical consistency of line identifica- 
tions, we modeled the negative Gaussian lines in the above fits with a photoionization model 
based on the code PHASE (Krongold et al. 2003). The PHASE code models self-consistently 
more than 3000 X-ray bound-bound and bound-free transitions imprinted by photoionized 
absorbers, given the ionization state of the absorber, its column density, and its internal 
turbulent motion. The parameters of the code are: 1) the ionization parameter of the ab- 
sorber U = - , where Q(H) is the rate of H ionizing photons, r is the distance to the 
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source, Uh is the hydrogen number density and c is the speed of light; 2) the equivalent 
hydrogen column density N H ; 3) the outflow velocity of the absorbing material V out , and 
4) the micro- turbulent velocity Vp ,. r h of the material. The abundances have been set at the 



Solar values (IGrevesse et al.lll993l ) 



A PHASE component (which we call system 1) with ionization parameter of logU = 
—0.60 ± 0.38 and column of logNH = 19.8 ± 0.2 cm 2 successfully reproduced the O VII Ka 
and K/3 lines at ~ 19.805 and ~ 17.08 A, at the blueshift of 32365 ± 335 km s" 1 with 
respect to source (Fig. 4). The addition of PHASE component "System 1" to Model-A 
significantly improves the fit (A% 2 = 15, A d.o.f. = 4). According to F-test, the absorber 
is present at confidence level of 99.87%. A second low ionization parameter component (we 
call this system 2) with logU = —1.2 ±0.21 is also required to fit other high velocity outflow 
absorption features such as of O VII and O VI , with outflow velocity of 31735 ± 275 km s^ 1 
(Fig. 4). Inclusion of this component improves the fit (A% 2 = 31, A d.o.f. = 4) over the 
previous model at a confidence level of more than 99.999%, according to F-test. The best 
fit PHASE parameters are reported in Table 2 and the best fit model are shown in Fig. 
5. Successfully modeling the residuals with two PHASE components robustly confirms the 
presence of high velocity outflows in Ark 564. 



5. Discussion 



5.1. Comparison with Models 



Several models suggest radiation pressure as the drivi ng mechanism to produce the 



typical low velocity outflows observed in Sey fert galaxies (jProga &: Kallmanl |2002| ; 12004 



pi 

2> 



Krolik fc Krisd 1 19951 . iDorodnitsyn et al.ll2008l ). But in radiation driven disk- wind models 
outflow velocities depend on AGN luminosity and th ese models cannot account for relativistic 
velocities in AGNs with Seyfert-type luminosities (IBarai et al.ll201ll ). Relativistic outflows 
in the UV have been d etect ed only in most l umino us broad absorption line quasars. Indeed, 
Laor &: Brandtl (120021 ) and iGanguly et al.l (120081 ) have shown that the maximum outflow 
velocity is proportional to L os , close to what is expected from radiation pressure driven 
winds. Observations of relativistic outflows in X-rays in moderate luminosity AGNs therefore 
pose intriguing puzzles; as shown in figure 6, the relativistic outflows (from iTombesi et al. 
20101 ) lie above the Ganguly et al. line of maximum velocity. The relativistic outflow we 
find in Ark 564, shown as a * in figure 6, also lies above the line. Thus these high velocity 
outflows appear not to be driven by radiation line pressure mechanism. 
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Kingl ( I2010l ) shock wind models produce winds with velocities v ~ 0.1c, but in quasars 
accreting at Eddington limits. In this model a high velocity ionized outflow collides with the 
ISM of the host galaxy, losing much of its energy by efficient cooling resulting in a strongly 
shocked gas. 

In the multi-dimensional Monte Carlo simulations of AGN disk-wind models of Sim et 
al., ou tflow velocities are as sumed to be at escape velocities, and not from ab-initio calcula- 



tions flSim et al.ll2008 



201 0J). It is pos sible that magnetic fields are important for launching 



such winds, as in jets ( jSim et al.ll2010l ); observations of relativistic outflows in moderate lu- 



minosity AGNs provide valuable constraints to the theory of magneto-hydrodynamic winds. 



5.2. Mass and Energy Outflow Rates 



In Paper-I we have showed that the typical warm-absorbers in Ark 564 (v out ~ 100 km s -1 ^ 
are not energetic enough to provide the effective feedback as required by theoretical models. 
Here we estimate the total mass and kinetic energy outflow rate of relativistic outflows in 
Ark 564 to determine if these can possibly be important for feedback. To measure the mass 
and energy outflow rates, we mus t know the location of the absorber which is measured only 
in a few sources (e.g. NGC 4051, iKrongold et al.l 120071 ); we do not know the location of the 
high velocity outflow in Ark 564 either. 

It is often assumed that the observed outflow velocity is the escape velocity at the 
launch radius r: i.e. r 



2GM BH 



"out 



There is no justification for this assumption, as shown 
by iMathur fc Stolll (j2009l ). No netheless, making thi s assumption provides us with a lower 
limit on the absorber location. iRomano et al.l ( 120041 ) determine the central black hole mass 



of Ark 564 to be Mbh < 8 x 1O 6 M . This leads to the minimum distance of the absorber 



of r r 



6.4 x 10 5 pc. For a sherically symmetric outflow, the mass outflow rate is M out 



Cf M v out (Ar) -1 where M is the mass in the outflow in a shell of thickness Ar and Cf is 
the covering fraction. Therefore M out = ^nCf r 2 Nh m p v out (Ar) -1 . Substituting Ar with 
r min and Cf — 1 gives us an upper limit on the mass outflow rate which is M out {r mir ^ < 
1.4 x 10 21 g s -1 . The corresponding kinetic energy outflow rate is E K = \M out v 2 ut , so we 
again have an upper limit Ek < 4.5 x 10 39 (r/r m j n )erg s _1 . This is a tiny fraction of the 
bolometric luminosity of Ark 564, approximately Ek/L^i <~ 0.001%. This upper limit is 
much sm aller than the ratio of E K /L hn i ~ 5% expect e d if the outflow is to b e important for 



feedback flSilk fc Reeslll998l . IScannapieco fc Ohll2004! . Idi Matteo et al.ll2005f ) 



The iKingi (120101 ) model for the quasi spherical shock wind produces a velocity gradient 
in which higher ionization elements are present at higher velocities. In such a scenario it is 
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possible that the column density we measure is not the column of the entire outflow, but 
only that of the shells producing the observed absorption lines; the total column could be 
much higher. We can estimate the mass loss rate of such a outflow by using the expression 
from King et al (2010): for a solid angle 4irb, mass loss rate is M out = 4nbnr 2 v out m p , where 
n and m p are the particle density and proton mass respectively. Mass conservation in the 
flow requires the product nr 2 v to be constant. As nr 2 = 4 ^* , the c onservation of mass 
in a r adial outflow yields a linear correlation of velocity and ionization ( jPounds fe Vaughan 



20111 ). For the range of b ~ 0.3 — 1, the mass outflow rate and kinetic power of the outflow 



varies from (1.3 — 4.3) x 10 25 g s _1 and (0.6 — 1.9) x 10 erg s _1 respectively. These estimates 
are close to the bolometeric luminosity of the source (18%-60%), much larger as required by 
feedback models. 



5.3. Comparison with other high velocity outflows 



Before this work, the high velocity outflows were detected mostly with Fexxv and/or 
Fe xxvi absorption lin es in hard X-ra y band (§1). In handf ul of quasar such as PG11211+143 
dPounds et al.ll2003ah . PG0844+349 (bounds et al.ll2003bh and MR 2251-178 JGibson et al. 
20051 ) high velocity outflows were detected at soft X-ray energies, but these detections were 
either based on absorption lines of low statistical significance or have strong contamination 
from the halo of our Galaxy. 

Here we present the detection of high velocity outflows in Ark 564; this is the first time 
such outflows are found in the soft X-ray band in a Seyfert galaxy. We firmly establish 
the presence of high velocity WAs first through identifying number of ionic transitions at 
similar velocity and furt her successfully modeling these features with photoionization models. 
Papadakis et al.l ( 120071 ) also detected an absorption line at 8.1 keV in the low resolution 
CCD spectra of Ark 564 and assuming that this line corresponds to Fe xxvi , they suggested 
the presence of highly ionized, absorbing material of logNu > 23 cm -2 outflowing with 
relativistic velocity of 0.17c. If the presence of such a feature and its identification are correct, 
then this is suggestive of a velocity gradient with higher charge states such as Fexxvi at 
higher velocity. Interestingly, the two "other" lines we identified in Ark 564 (§3.0.3) also 
have somewhat higher velocity than System 1 and System 2 and have higher charge states. 
This is exactly as expected from the models of iKingl (120101 ) . 
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6. Conclusion 

We report on the discovery of high velocity outflows in the NLS1 Galaxy Ark 564. These 
absorbers are identified through multiple absorption lines of Ovi , Ovil Ka, Ovu K/3, 
Ovill and Fexvil at blueshift of ~ 0.1c (with respect to the source) detected in the 
Chandra HETG-MEG spectra. The two observed velocity components are well fitted with 
two photoionization model components. Both absorbers have low ionization parameter of 
log£ = -0.60 ± 0.38 and -1.2 ± 0.21 and low column densities of \ogN H = 19.8 ± 0.2 and 
20.0 ± 0.3 cur 2 and are required at high significance of 99.87% and > 99.99% respectively. 

Without knowing the location of the absorber, its mass and energy outflow rates cannot 
be well constrained; we find E(outflow)/L bo ranging from < 0.001% to 60% using different 
assumptions. Determining the absorber location is therefore very important for providing 
meaningful constraints. This can be achieved through studying the response of absorption 
lines to continuum variations. This is the first time that absorption lines with ultra-high 
velocities are unambiguously detected in the soft X-ray band. The presence of outflows 
with relativistic velocities in AGNs with Seyfert-type luminosities is hard to understand and 
provides valuable constraints to models of AGN outflows. Finding such relativistic outflows 
in several other AGNs and measuring their mass/energy outflow is therefore important. 
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Table 1. Absorption lines identified with relativistic outflows in the HETGS Spectrum of 

Ark 564. 



Ion 


^obs 




Vout 


EW 




A 


A 


km s _1 


mA 


Ovi 


20.25 


22.05 


30964 ± 355 


16.1 ±3.4 


O vii Ka 


19.84 


21.60 


31042 ± 362 


16.5 ±2.4 


O vii Ka 


19.80 


21.60 


31588 ± 362 


15.6 ±2.5 


O vm 


17.35 


18.97 


32112 ±410 


11.8 ±2.7 


O vii K/3 


17.08 


18.63 


31634 ± 416 


4.0 ±2.4 


Fe xvii 


13.62 


15.02 


34951 ± 316 


5.0 ± 1.7 



Table 2: Model parameters for Relativistic Outflows 



Parameter System 1 System 2 

Log U -0.60 ±0.38 -1.2 ±0.21 

Lo#A^(cm- 2 ) 19.8 ±0.2 20.0 ± 0.3 

V Turb (km s" 1 ) 85 ± 14 92 ± 12 

Vout (km s" 1 ) 32365 ± 335 31735 ± 275 

Ax 2a 15 31 
a Improvement in x 2 to fit after adding model component 
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Fig. 1. — Residuals of data: Model- A fit to the coadded spectrum of Ark 564, in the ob- 
server frame (Model-A: Two component WA model plus continuum, emission lines and local 
absorption, from Paper-I). The possible transitions due to high velocity WAs are indicated 
with red labels. 
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Fig. 2. — Residuals of data:Model-A fit to individual spectra from the 2008 observations. 
The red boxes mark the regions near 19.804, 19.844 and 20.24A where the strongest lines 
are detected in the combined spectrum. We see that the same lines are detected in the 
individual spectra as well. 
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Fig. 3. — The Ark 564 MEG spectrum fitted with Model-A plus six Gaussian lines. The 
line identifications of high-velocity outflow components are indicated with blue labels. The 
spectrum is presented in the observer frame. 
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Fig. 4. — The high velocity outflow absorption features fitted with two PHASE components 
of outflow velocity of 0.103c (red) and 0.105c (blue). The ionic transitions are labeled in 
black. The spectrum is presented in the observer frame. 
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Fig. 5. — The two component PHASE model, one with v out = 0.105c (blue) and other with 
v ut — 0.103c (red), showing the enormous blueshift of absorption lines with respect of source 
redshift. 
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Fig. 6. — Outflow velocity plotted as a function of AGN l uminosity. The solid line repre- 
sents the upper envelope relation from lGangufy et all ( 120081 ) modified to plot the bolometric 
luminosity instead of the 3000A luminosity. The points are for the relativistic outflows in 
Tombesi et all (l2010f ). Ark 564 ultra-fast outflow is shown by a star (this work). It is clear 
that these ultra-fast outflows are not confined within the Ganguly et al. envelope. 



